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The non-equilibrium time evolution of an Anderson quantum dot is investigated. The quantum dot is coupled 
between two leads forming a chemical-potential gradient. We use Kadanoff-Baym dynamic equations within 
a non-perturbative resummation of the s-channel bubble chains. The effect of the resummation leads to the 
introduction of a frequency-dependent 4-point vertex. The tunneling to the leads is taken into account exactly. 
The method allows the determination of the transient as well as stationary transport through the quantum dot, 
and results are compared with different schemes discussed in the literature (fRG, ISPI, tDMRG and QMC). 
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I. INTRODUCTION 

Due to the recent advances in nanoscale technology, it has 
become possible to study electron transport experimentally in 
nano-devices such as single molecules, artificially designed 
quantum dots or nanotubes [[ll-lll]. This naturally incited new 
interest, and quantum dot models have been the subject of 
much theoretical effort recently. The quantum dot also serves 
as a good playground for theoretical methods: it allows for 
many-body effects as well as non-equilibrium phenomena in 
a regime where linear response theory no longer applies. The 
description of co-tunneling processes such as the Kondo ef- 
fect and the emergence of the Kondo scale have proven to be 
especially challenging. This is the case, in particular, for the 
exponential coupling dependence of the Kondo scale. 

The Kondo regime has been investigated perturbatively us- 
ing Fermi liquid theory exploiting integrability of the An- 
derson model 101, or the non-crossing approximation at in- 
finite local coupling U Ht). Several renormalization-group 
methods have been used to investigate the stationary state of 
the system: the perturbative real-time renormalization group 
(RG) J^, non-equilibrium extensions of perturbative RG lIlOl - 
\VM . and the functional RG approa ch in its generalization to 
non-equilibrium situations |[l3l - [l9ll . Bethe-Salpeter equations 
in so-called parquet approximation ll20ll were analysed with 
respect to the Kondo resonance in Refs. ll2m23ll . The numer- 
ical renormalization group (NRG) method has also been suc- 
cessful in describing such impurity systems, and it has been 
generalized to describe time-dependent non-equilibrium sys- 
tems iiHi. 

Another line of thought is represented by the quantum 
Monte Carlo method (QMC) ll28l] . which is numerically exact, 
but only short simulation times are possible since the method 
is plagued by the sign problem at small temperatures. 

The transient behaviour of the system as well as the station- 
ary state have also been explored by methods such as the time 
dependent density matrix renormalisation group (tDMRG) 
II29I - U2II which allows the simulation of the time evolution of 
pure states, and the iterative real-time summation of the path 
integral (ISPI) |[33| which is numerically exact, but depends 
on the correlation time of the system being small. Several 
of these theoretical methods were recently compared in ll34ll . 
which also contains a concise list of related previous studies. 



In this article we use the two-particle-irreducible (2PI) ef- 
fective action or O-functional ||35| - |37I1 . to derive Kadanoff- 
Baym equations of motion describing the transient and sta- 
tionary transport through a quantum dot. The power of this 
method is that symmetries of the system are conserved during 
the time evolution, i.e., particle number and energy conserva- 
tion is satisfied. This method has been successfully applied to 
the thermalisation of relativistic and nonrelativistic systems of 
bosonic and fermionic gases Issl - liill . 

We apply the 2PI method to the single-impurity Anderson 
model of a quantum dot within two different truncations: the 
loop expansion up to two-loop order (in the self-energy), and 
a resummation, which includes the contribution of the spin- 
aligned bubble chains summed to all orders. As shown below, 
the effect of this s-channel resummation can be seen as using 
a frequency dependent 4-point vertex. The leads can be taken 
into account exactly in this formaUsm: integrating them out 
they give a contribution to the self-energy of the dot-electrons. 

After benchmarking our method against exact solutions for 
an isolated quantum dot we compare our results for the tran- 
sient and stationary behavior of the dot coupled to the leads 
with results found in other schemes, including fRG, ISPI, 
tDMRG, and QMC calculations. While finding good agree- 
ment in the long-time stationary limit with the behavior found 
with these schemes our approach also allows to treat the tran- 
sient build-up of the current through the dot. Comparisons 
of the stationary current for different interactions show good 
agreement with the alternative schemes over a large range of 
interaction strengths. These results underline the necessity of 
a nonperturbative resummation. 

The structure of the paper is as follows: The Anderson 
model and its ingredients are defined in Section|II] In Section 
|lll]we summarize the formulation of the 2PI eff'ective-action 
approach for nonrelativistic Fermions on an isolated quantum 
dot. We introduce the formulation for the dot coupled to leads 
in Section |IV] In Section |V] we compare our method with 
exact results, for cases where an analytical treatment is possi- 
ble. In Section |VT] we present our results for the build-up of 
the current and for its stationary characteristics, and compare 
these to results in the literature. We draw our conclusions in 
Section |Vn] 
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Figure 1: The schematic representation of the energy levels in the 
Anderson dot system, for the zero magnetic field case. The lower line 
of the dot represents the energy level of the first electron occupying 
the dot, the upper line is the level of the second electron occupying 
the dot. 



II. THE ANDERSON MODEL 



The Anderson Hamiltonian describing the electrons on the 
quantum dot and on the leads reads 



H - //dot + //leads + //tunnel - ^ /^Oo-Wtr + Utljni + 

+ X ^f^P^pa^kpcr - Yj^tpcl^^J^ + fpdlckpa) ( 1 ) 
kpcT kpcr 



tonian ([T]i are given by 

5dot = J" dt^dl(idt-EQa-)da-Ud^^d^d^^di, (2) 

5lead.s = I dt^cl^^^iid, - ekp)Ckpa; (3) 
kpa 

■Stunnel - \ dt^ {tpc[ da- + fpdlckpa)- (4) 

Our most important observable is the time dependent cur- 
rent through the dot 

^^'^ " -? Z {P'p<''lp.''-^ - PC(dUkpa)) ■ (5) 

kpcr 

This can also be written as / = (Il-Ir)/!, where Ip - -eNp{t) 

and Np{f) — {Yjkcr '^pjto-'-p^o"^ '■^^ number of electrons on the 
leads. The stationary current can simply be obtained by wait- 
ing for the transient behavior to die out, such that the system 
is sufficiently close to the final, stationary state. 

An important scale of the physical processes described by 
the Anderson model is given by the Kondo temperature. This 
temperature marks the onset of the Kondo effect which is due 
to the formation of singlet states of intinerant and localised dot 
fermions and characterised by a rising resistivity of the dot at 
low temperatures. The Kondo temperature is found by simple 
renormalisation-group arguments to be 



r-^Vf -plai- 



ce) 



for the particle-hole symmetric system, where £0 - -U jl. 



where cr - +1/2 ^ f, X is the spin index, p = + ^ LjR labels 
the leads on the left and right, respectively, and k is the index 
of the spectrum of the lead electrons. The leads on the left 
and right have a chemical potential jip - peV/2 which models 
the bias voltage. For the determination of the energies etp, 
see Sect. |IV] The occupation-number operator on the dot is 
written as: «o- - dldg-- The dot electrons have an interaction 
term ~ U due to Coulomb repulsion. The one-body energy is 
£00- - E{) + crB. The first term is controlled experimentally 
through the gate voltage whereas the second (Zeeman) term 
by the magnetic field B. The energy levels of the system are 
represented in Fig.[T] 

The tunneling between the dot and the leads is quantified by 
the tunneling strength tp. We will assume that the leads are in 
thermal equilibrium at some temperature T and use the usual 
wide-band limit with a constant density of states p around the 
Fermi surface. In this work we use a symmetric tunneling rate 
T - ti - tR. The dimensionful quantities of the system can 
be expressed in units of the hybridization F — 2n\T\-p which 
quantifies the dressing of the dot by the leads. 

We will use, in the following, a Lagrangian or action for- 
mulation of the Anderson quantum field model. The contribu- 
tions to the action corresponding to those forming the Hamil- 



III. 2PI FORMALISM FOR THE ISOLATED DOT 

In this paper, we study Kadanoff-Baym-type dynamic equa- 
tions for correlation functions describing transport through the 
Anderson quantum dot. We derive these equations using the 
two-particle irreducible (2PI) effective-action formalism, also 
known as <I>-derivable approach. In the following, we sum- 
marize known basics about the 2PI effective-action approach 
ll35l - [37ll and introduce specific details of its implementation 
for the Anderson model. 

Since we will be interested in initial value problems, specif- 
ically, in the evolution of multi-time correlation functions 
starting from values given by some initial state, it is conve- 
nient to work in the Heisenberg picture. In the correspond- 
ing functional-integral formulation, the time integrations in 
the action are defined to run along a Schwinger-Keldysh con- 
tour C, which goes from the initial time fo to some final time 
t, and then back to fo i45u461. 

The KadanofF-Baym equations will yield the evolution of 
the time-ordered two-point function: 

D^At, t') = 0c(r - t'){d^{t)d\{t')) - ©c(f' - t){d\it')d^(t)) 

(7) 
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Figure 2: (color online) 2PI diagrams of the loop expansion ofTjlD]. 
(a) The two lowest-order diagrams of the loop expansion described in 
sections fill Al and IlII B I with each vertex coupling up spins to down 
spins. Black dots represent the bare vertex ~ U, blue (green) solid 
lines the up (down) spin propagator Do-o-, o" =1, 1- (b) Diagram rep- 
resenting the resummation approximation explained in Sect. IIII CI 
The wiggly line is the scalar propagator Gu . This propagator is rep- 
resented as a sum of bubble-chain diagrams, see OSt . The same 
diagram with down spins also appears in r2[ZJ]. 



where &c(t - t') is a ©-function on the Schwinger-Keldysh 
contour, and evaluates to 1 (0) if t is later (earlier) to t' . 

We begin with the simpler case of the isolated quantum dot. 
Without leads the 2PI effective action IBTIl reads 



r[D] = -/Tr[lnD ' + Dq'd] + r2[D] + const. 



(8) 



Here r2[D] can be written as the sum of all closed 2PI di- 
agrams constructed from bare vertices and full propagators. 
Two-particle irreducible are those diagrams which do not fall 
apart upon cutting two lines IstIi . As an example we show, in 
Fig 12^, the two diagrams of lowest order in the bare coupling 
U, represented as a black dot. Full lines denote different spin 
components of the full propagator or two-point function D. In 
(IS), the free inverse propagator is given by 



(9) 



The Kadanoff-Baym equations of motion are calculated from 
the stationarity conditions for the action 



6Y[D] 

6DAa(t'j) 



0. 



(10) 



This equation gives the well-known Dyson equation 
D^\(t,t') - DoJ^^(f, f') - So-,i(f, f'), where the self-energy is 
given by 



_6r2[D]_ 



(11) 



We use the decomposition of the time-ordered two-point func- 
tion ^ 



Da-Ait, t') = FcrAit, t') - -pa-A{t, t')signc(t - f), (12) 



where 



FaA(t,t') = ^{[dAt),dl(t')]), 

Pa-Ait, t') = i{{dAt),d\(t')}), 



(13) 



are the so-called statistical and spectral correlation functions 
in real-time representation. signp(f - f) - -I + 20c(f - t') is 
a sign function on the contour For the equal-time arguments 
the following identities hold 



FaAt,t)^^ 



■ nAt), Po-A{t, t) = iSaA 



(14) 



where «o-(f) is the mean number of fermions occupying the dot 
level with spin cr at time t and the equation for p(f, t) follows 
from the fermionic equal-time anticommutation relations. The 
symmetry relations 

F^Ait, t') = FaM, t)\ PAcAt, t') = -paA(t', tf. (15) 

are a direct consequence of the definitions of F and p. For 
its numerical solution it is convenient to rewrite the Dyson 
equation for D(t, f) in terms of two equations for F{t, f) and 
p(f, t') whereby the time integrals over the contour C are re- 
placed by simple integrations along the real-time axis. In this 
way one arrives at the usual 2PI or Kadanoff-Baym equations 
of motion for the complex functions F{t, t') and p{t, t'). 

We will assume throughout that correlations between up 
and down spins vanish, such that Dua ~ S,rA- We also in- 
troduce the notation: D^,r s (no summation!). With this, 
the equations of motion read 

(id, - Ma-it)) Pa-it, t') ^ J di^ ^ait, u)paiu, f), 

iid,-Mait))FAt,t') = r dulP,it,u)Fpiu,t') 
Jo 

duLl.it, u)pa-iu,t'), (16) 



Jo 



where the decomposition 

LAt,t') - -nt\t)5it-t') + 



(17) 



+E^(f,f')--signc(r-f')5:^(f,0. 



of the self-energy into local, F and p type terms has been used. 
The time-local energy term in ( fTSI l is given by 



M^it) = + l.f\t) 



(18) 



and includes the mean-field shift originating from the double- 
bubble diagram in Fig |2^. The above integro-differential 
equations are equivalent to the exact Kadanoff-Baym equa- 
tions and include higher-order correlations through the non- 
Markovian memory integrals on their right hand side. 

Including the second-order 'basket-ball« diagram in Fig|2^ 
beyond the mean-field approximation yields a perturbative 
small-coupling approximation of r2[Z)], 



Y2[D]^-U ^^dtDAt,t)DAt,t) 



+ ^U- £dtdt'DAt,t')DAt',t)Diit,t')Diit',t). (19) 

In the self-energy this corresponds to taking into account the 
tadpole and sunset diagrams. 
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A. Mean-field approximation 

The one-loop contribution to the self-energy gives the 
■field approximation of the dynamic equations ( fTSI l, 

(20) 



mean 



YJ^\t,t')^iU6{t-t')D^At,t). 



Note that D,j-(f, t') has a jump at equal time arguments. By 
analyzing the path-integral construction one finds that using 
the specific operator ordering in the Hamiltonian (IA2I) implies 
that the self-energy takes the form 



(21) 



See AppendixlAlfor details. The resulting equations of motion 
to this order of approximation are local in time and therefore 
Markovian. They constitute the 1-loop or mean-field approx- 
imation. 



B. Second-order approximation 

We finally derive the order- 1/^ contribution to the self ener- 
gies, also called second-order Born approximation, 



i^Sunset 



(22) 



or, in components. 



pa(\F-j-^\p-jy 

+Fa{F-^p*_^ + F*_^p-a-)\ 



(23) 



where we have omitted arguments. The contributions (l2Tl i 
and ( |23] | are summed up to be used in Eq. ( fT6] l. with (l2Tl i 
contributing to S*^*. Some details concerning the numerical 
implementation of the above equations are given in App. |B] 



C. s-channel resummation 

In this section we will summarize the i-channel resum- 
mation which leads substantially beyond the coupling ap- 
proximations discussed before. Specifically, this involves a 
summation of the bubble chain contributions with alternat- 
ing spins where, in each bubble, the two propagators de- 
scribe the same spin component. This is similar to the next- 
to-leading-order 1 IN approximation for A^-component scalar 
fields llssl I49I1 . An elegant way to perform the resummation 
involves a Hubbard-Stratonovich (HS) transformation. Since 
the interaction vertex couples 'up' spins with 'down' spins, 
the bubbles in the chain have alternating spins. The action 

■5 dot = J" ^/f ^ dliid, - Eoa-)da - Ud!^d^d'^di (24) 



is rewritten using auxiliary scalar fields xi and xi by use of 
the substitution 



J A- 'J -^xtAx + 2JtX 



(25) 



where 



1 

'-2 



( d^^df \ 



d[di 

The resulting action reads 

SAoAdcr,dl.,Xi\ 



A = 



1 /O 1 



2f/ \ 1 



(id, - Eoa-)da- 



+JJX1X2 + d}^d-^xi + dldixi- (26) 

The free inverse propagators are read off from the quadratic 
part of the action 

iGQ\t,t') = 2A6{t-t'), 

iDQl{t,t') = {id,-Eo^+Xi)S{t-t') (27) 

where the free propagator Go of the scalar fields is a 2 x 2 
matrix. Accordingly, we call G the propagator of the scalars, 
and Xi is the one-point function or expectation value of the 
auxiliary fields. 

The corresponding 2PI effective action can be written as: 

T[G,D,x] = SUdl = flf^ = Q,x] - !Tr[lnD-' + 

H-^Tr [inG ' + Gg 'g] + YiiD, G] + const. 

(28) 

where T2[D,G'\ contains all closed 2PI diagrams built from 
the 3-point vertices of the action ( l26t and full scalar and 
fermion propagators. The lowest-order contribution is shown 
in Fig.|2]3. 

The stationarity conditions give the Schwinger-Dyson 
equations. In this case, the field average of the scalar fields, in 
contrast to that of the fermionic fields, is nonzero, so we have 
stationary equations of the form 6TI6x - 0- The resulting 
equations read 

Xxit) = f/Du(f, f), Ut) = UD^^iu t), (29) 
D^^D = 2 * D + 5, (30) 



G-^G^n*G + 6 



(31) 



where we again suppressed the time arguments, and * stands 
for convolution on the contour C, 



(A * B){t, f 



') = j dzA{t,z)B{z,t'). 



(32) 



S and n are the self-energy of the fermion and the boson 
fields, respectively: 



n(f,f') = 2/ 



:dr2[D,G] 
dDAt'J) 
dr2[D,G] 



-DAt,t')G^At,t'), 



dG{t,t') ' 

= DAt,t')DAt',t), 



ni2 = 021 = 



(33) 
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where 1-1, J, = 2 is used for the field indices. 

The equations concerning the scalar fields are constraint 
equations, which do not contain any time derivatives, because 
thex fields are auxiliary, non-dynamical fields. The r2 part of 
the action to lowest order is obtained as (See Fig. |2j3) 



T2[D,G] 



2/ 



y)D^(y,x)G^Ax,y). (34) 



From the constraint equation (IJTT i for G, one can see that 



(35) 



- iU 
+ U 



n22*nii 



nii*n22 
4/ 1122 * rill * 1122 







iiii * n22 * 111 



where we have again omitted the (f, f') arguments. Inserting 
Gil and G22 into the self-energy Z of the fermions in eq. ( [33T l, 
one finds a sum of bubble chains with alternating spins being 
generated as mentioned initially. The decomposition of the 
propagator G and the self energies into statistical and spectral 
parts is deferred to AppendixICl 



IV. LEADS IN THE 2PI FORMALISM 

In this section we couple the single dot system to the leads 
taking into account the terms Q and (HJi in the action. We can 
calculate the leads' contribution exactly, since ([3j and (|4]i only 
include terms up to quadratic order, by integrating these out. 
The initial density matrix of the full system is assumed to be 
of the product form 



p(t -0) - pdot ® Pleads, 



(36) 



i.e., no correlations between dot and lead electrons exist ini- 
tially. Pdot is described in Sect. IV Al while pieads describes 
a grand-canonical ensemble, as shown below. To include a 
thermal density matrix for the lead electrons, we expand their 
time-contour with a vertical part after C from fo to to - i/3 for 
each momentum mode of the lead electrons. Then, the contri- 
bution to the dot degrees of freedom is 



,5env = -/|t|2 J dt dt'd\t)S-\t, f) d(t') (37) 



where S (t, t') is the free action of the lead electrons, 

S(t,t')^{id,-e(t))6{t-t'). (38) 

Note that we have chosen symmetric tunneling t - ti - Ir. 
Here f is a complex variable, which lives on the extended 
Schwinger-Keldysh contour. e(t) = ekp if t is real, with k and p 
being the indices of the mode considered, and e(f) = ekp- 1^ if 
t is on the vertical part of the contour This prescription gives 
the correct initial density matrix corresponding to the grand 
canonical ensemble. 



To obtain the contribution to S^m, we have to solve the 
equation 

{id,-e{t))A{t,t')^6c{Ut'), (39) 

using the antiperiodic boundary condition A(0, t') - 
-A{-i/3, t'). The result can be written in the form 

S-\t,t') = A(f,f') 

= A>(f, f')0c(f - t') - A<(f, f')0c(f' - t), (40) 

A<(f,r') = -if{ekp-^l)e-•'^"'■'-''^ (41) 

where f{x) - 1/[1 + exp(j6x)] is the Fermi function. After 
integrating the lead electrons out, this appears as part of the 
free propagator of the dot electrons, so we have 



[Do' + /a] * D = S * D + 5. 



(42) 



After decomposition one finds that the effect of one lead- 
electron mode can be seen as a contribution to the self-energy 
of the dot electrons if and If: 



(43) 
(44) 



One can also argue the following way: If the tunneling terms 
are used as a vertex, they give a one-loop diagram in r2 [£>] , 
half of the loop is the propagator of the electron on the dot, 
the other half is the propagator of the lead-electron. After 
taking the derivative of V{D] one sees that up to a factor the 
contribution to the self-energy is the propagator of the free 
lead electron. This gives the same result as above. 

Hence, the contribution only depends on the difference of 
the time arguments. This is to be expected, since the lead 
electrons are in thermal equilibrium. Now we have the contri- 
bution of one lead particle. The wide-band limit corresponds 
to integrating over a continuum 



^lead - 

J-D 

with p(e) - const. With ( l44l i this leads to 



If 

lead 



f 

J-D 



de iM^e-'" 



-ilrn- 



sin Dt 



t 



It is useful to consider the limit £) — » 00, where 

^Ld(f'0 = -2'^M''^(f-0- 



(45) 



(46) 



(47) 



The calculation of 1.'^ can only be done analytically at zero 
temperature. The finite-temperature correction is calculated 
numerically, whereas the zero-temperature contribution is a 
principal value. 



L,R 



-iftt 



(48) 
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To calculate the time dependent non-equilibrium current, 
one adds a further source term to the action: 

^ " = ? Z P^'p'lp.'^- - f^dUkpaXt,,,) (49) 

Then the current at some measurement time t,„ can be written 

as 



d 

I{t„,) = -i— In Z[7j] 
orj 



(50) 



,,=0 



After integrating out the leads, one gets the following contri- 
bution to the action: 



iS lead - -i 



Y^JdtJ dt'd'(t)I^Ap(t, f) 



+ '-^Apit, t') [-6{t - t,„) + 6{t' - f,„)] \d{t') (51) 



where we dealt with the first term in the previous subsection, 
see Eq. ( l40b . The second term will provide a contribution to 
the current which thus reads 

hit) = e^i r dt[Al(t„,t)pAtJm) 

- A^(f„„ f)F^{t, f„,) + Af(f„„ t)pl{t, tj 
-Al\t„„t)F*^(t,t,„)]. (52) 

For the full current one needs I - {I^ - Ir)/2. If the leads are 
symmetric (tunneling rate, temperature and level density are 
the same for the left and the right lead) the Ap terms cancel, 



'symm civv^ 



dt[A'[(t,„,t)-A^{t„,t)\pAt,t,r,)- (53) 



V. COMPARISON WITH EXACT RESULTS 
A. Exact time evolution without leads 

Without the leads, the time evolution of the dot is a sim- 
ple quantum-mechanics problem of solving the Schrodinger 
equation. To compare with results obtained with the above 
functional methods, we need to choose an initial density ma- 
trix for which the 1 -point, 3-point and higher connected n- 
point functions are zero: Tiipda-) — 0, Tiipda-dl^d-a-) - 
0, Tr{pF,jA) = So-aH/^ - rig-), TT(pdo-d-o-) = 0, etc. These 
requirements fully constrain the density matrix. 

Then one can solve the time evolution analytically by diag- 
onalizing //dot defined in Eq. ([T]). The commutator part of the 
2-point function is given by 

Fn(0,t) = Fl^(t,0) = ^ [(1 -«i)(l -«2)e'^»>' 



2 

-H (1 - ni)n2e'^^°'^'^^' - «i(l - nz)^'^"'' - «i«2e'<^°'^^*' 



]■ 

(54) 
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Figure 3: Time evolution of the two-time congelation function 
Fn(0, for an isolated quantum dot: Comparison of 2PI and exact 
results obtained by directly solving the Schrodinger equation. The 
magnetic field is zero. While the exact evolution shows revivals, the 
oscillations are damped out within the 2PI approach. 
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Figure 4: The same time evolution as in Fig. [3] Comparison of 2PI 
and exact results. The difference between the exact and 2PI results, 
as well as the difference between the Taylor-expanded exact and the 
exact results are shown. 



We compare this non-equal-time two-point function to the so- 
lution of the Kadanoff-Baym equations obtained from the 2PI 
effective action. As one sees in Fig. [3j the 2PI solution is in 
good agreement in the begin ning , but does not show revivals 
as expected, cf., e.g., Refs. ISOl l5lll. Coupling to the leads 
causes a real damping which the 2PI approach is capable to 
describe. In Fig. |4] we plot the difference between the exact 
and the 2PI results for Fn{0, t), as well as the difference be- 
tween the exact and the Taylor-expanded exact results. This 
illustrates that the 2PI equations yield a well-behaved dynam- 
ics. In contrast to this, perturbation theory shows secularity 
problems. 
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Figure 5: The stationary cureent l l55t tiiroughi tlie quantum dot at 
U = 0, as a function of the bias voltage. The comparison of the 2PI 
results and the results of Ref. IT4I1 . which in this case are both exact, 
serves as a benchmark test. 
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Figure 6: Comparison of the transient currents obtained with the 
2PI equations and with the tDMRG method, at the symmetric point: 
Eo = -UI2 . The parameters are: {/ = SF.T = O.IF. Note that 
the different initial slopes are due to differently chosen initial condi- 
tions. While, in the tDMRG calculation, initial correlations between 
dot and leads are present, there are no such correlations in the 2PI 
approach. 



also give hints to the limitations of the .^-channel resumma- 
tion. We consider first the transient evolution of the current, 
showing the build-up of the stationary state and then analyse 
the resulting asymptotic current. 



B. The current atU = Q 

Another quantity which can be calculated exactly is the cur- 
rent at vanishing on-dot coupling U. This is due to the fact that 
at zero coupling the dot degrees of freedom are also quadratic. 
For details of the calculation, see, e.g., Ref. H. The result is 

I^- { de{f(e-^iL)-f(e-^iR))- (55) 

where f{x) = 1/[1 -1- e\p(J3x)] is the Fermi function. 

In Fig. |5] the currents as obtained from the analytical for- 
mula ( fSSl ) and from the 2PI calculation (after waiting for the 
stationary state), for U - 0, are compared. The agreement 
benchmarks our calculation, because we take the leads into 
account exactly, and atU = the 2P1 treatment is also exact. 

VI. TRANSIENT AND STATIONARY DYNAMICS 

In the remainder of this paper we present our results for the 
transient build-up of the current through the quantum dot, for 
different coupling strengths, bias voltages, and temperatures. 
We compare our results with those obtained with alternative 
methods as given in the literature. As we will show, these un- 
derline the necessity to go beyond the perturbative coupling 
expansion of the 2PI part r2 of the effective action. Our results 



A. Numerical comparisons of different approaches 

For the comparisons we use data that was published in 
Refs. iHH. For details about the fRG (functional Renor- 
malisation Group), ISPI (Iterative Summation of real-time 
Path Integrals), QMC (Quantum Monte Carlo) and tDMRG 
(time-dependent Density Matrix Renormalization Group), 
see, e.g., Ref. ll34ll . and references therein. If not otherwise 
stated, we will use the particle-hole symmetric point, where 
Eo = -U/2, and U ^2r,B ^0,T ^ O.lf. 

In Fig. |6] we compare the time dependence of the current 
through the dot as obtained with the full 2PI equations of mo- 
tion with the results of the tDMRG method. The stationary 
values obtained with the two methods are close to each other, 
but the time dependences are different. The initial slope of 
the cuiTent for the 2PI curve is given by Ih/{e^V) - 2Vt, as 
can be seen by evaluating the Heisenberg equation of motion 
/ = i/fi[H, I], with the uncorrelated density matrix (l36T l at the 
initial time. The tDMRG method uses an initial density matrix 
containing correlations between dot and lead electrons, which 
explains the different slope at the initial time. 

In Fig. |7]we compare, for small couplings, the stationary 
cuiTent with the corresponding results of the ISPI approach 
discussed in Ref. ifssll . The perturbative results were presented 
in Ref. |52i|. Note that we plot the interaction correction, that 
is, the difference of the currents for zero and non-zero cou- 
pling which, for the chosen couplings, is a small quantity. The 
Coulomb-blockade physics suggests that the sign of the cor- 
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Figure 7: The modification of the asymptotic stationary current due 
to nonvanishing interactions is shown as functions of the bias voltage. 
The results obtained with the 2PI method are compared with those of 
the ISPI calculations given in Ref. (33ll . and with perturbative results 
fromRef. (sl. 



Figure 8: The time dependent current through the dot is shown for 
different couplings U, bias voltage V, and temperature T. We com- 
pare results obtained within the second-order truncation of the 2PI 
effective action with those derived after j-channel resummation. For 
U = 2T there is perfect agreement. Note that, for the bigger cou- 
pling, the Kondo temperature ^ is Tk = 0.33r such that the system 
is already in the Kondo regime. 



rection should be negative, as transport is suppressed for non- 
degenerate levels on the dot. One observes that the agreement 
between the results of the two methods, within the range of 
couplings considered, is very good. For U = O.IF the pertur- 
bative results are close to the ISPI and 2PI results. We find 
that for U - 0.3r the perturbation theory is of limited accu- 
racy while the 2PI and ISPI results agree well. 

In Fig. [8] the time dependence of the current is shown for 
different choices of the interaction, the bias voltage, and the 
temperature. We compare the second-order coupling expan- 
sion, see Sect. IIIIBI with the ^-channel resummation intro- 
duced in Sect. IIIICI For the case of stronger interactions 
we compare the stationary current reached with the result ob- 
tained with perturbative RG methods, as given in Ref. ll33ll . 
While for small coupling the truncations agree and reproduce 
the ISPI result given in Ref. llssll . deviations appear for larger 
U, where the system is expected to enter the Kondo regime. 

In Fig.|9]we compare the stationary current for bigger cou- 
plings, at the particle-hole symmetric point [Eq - -U/2), with 
corresponding results from the other methods. With this we 
provide an add-on to the comparison presented in Ref. ll34ll . 
One finds that for U = 2r all methods agree. At f/ = 4r 
deviations start to appear which become larger at U - SF. 

In Fig. [TO]we test the mixed-valence regime, by comparing 
the steady-state current, using Eq - F = -U I2 + 2Y. For small 
bias voltages one sees good agreement between all methods. 
For further details concerning the other methods see the dis- 
cussion in Ref. [[3lll . We furthermore compare our results with 
those from the fRG and ISPI methods at non-zero magnetic 
field, as shown in Fig. [TT] and find good agreement for the 
full range of bias voltages. In Fig. [12] we finally compare the 
temperature dependence of the steady-state current with fRG 
and ISPI results, for U ^ 2Y, T = O.IF. One finds good 
agreement for all three methods, with the largest deviations in 
the high-temperature large-bias regime. 

Overall, we conclude that in the moderate-coupling re- 
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2PI 
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Figure 9: Comparison of the results obtained with different methods, 
at the symmetric point Eq = -U/2. The stationary current is shown 
as a function of the bias voltage for several coupling strengths. The 
2PI and the ISPI methods use T = O.IF. 



gion, the 2PI method gives reliable results irrespective of the 
bias voltage and the temperature, both at the symmetric point 
(£0 = -U/2) and in the mixed- valence regime. 



B. Effective coupling for ^-channel resummation 

We close this section with an analysis of the effect of the 
i-channel resummation on the effective coupling strength on 
the dot. This resummation allows us to write the equation for 
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Figure 10: The stationary current in the mixed-valence regime for 
the coupling of = 2r , and Eg = F. See Ref. jsill of a detailed 
discussion of the other methods. 



Figure 12: The stationary current is plotted as a function of the bias 
voltage for several temperatures. 
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Figure 11: The stationary current as a function of the bias voltage, 
for several magnetic field strengths. 
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Figure 13: The effective coupling for various couplings as a function 
of the frequency. 



the fermion self-energy as 

1.'" = U(FA'' -^pAP)^UU,«{Fn'' -^pRP), 



where 



1 -P 



R|2 



l(n«)2-i|: 



(56) 



(57) 



is a real, w-dependent effective coupling function. Details of 
its derivation are given in AppendixlD] 

Note that, before resummation, if and If had the same 
form, with U^s replaced by U, see in Eqs. ( l23T l and ( IC2b . The 
time arguments in ( l56l ) are always (f, f). 

One can thus understand the effect of the resummation as 
introducing a frequency dependent 4-point vertex. This vertex 
only depends on one frequency, because only the i-channel di- 
agrams are resummed. In Fig. [13] we show the effective cou- 



pling for various bare couplings U. For small couplings, one 
has J/eff ~ 1 because 11 is proportional to the coupling U . 
For medium-sized U, the effective coupling shows big devi- 
ations from 1, and for large U the effective coupling again 
approaches one. 



VII. CONCLUSIONS 

In this paper we have studied the real-time evolution of 
the Anderson quantum dot, using the two-particle-irreducible 
(2PI) effective-action approach. The dynamical evolution of 
correlation functions is described by Kadanoff-Baym equa- 
tions of motion derived from the effective action in nonper- 
turbative approximation. Two different approximations were 
considered: the second-order (sunset) coupling-expansion ap- 
proximation and the ^-channel resummation approximation, 
in which the bubble chains (bubbles with equal-spin propaga- 
tors, alternating spins in adjacent bubbles) were summed to 
all orders. 
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The leads were taken into account by integrating them 
out exactly, using a grand-canonical initial density matrix. 
At nonzero bare interaction strength U of the dot fermions 
we compared our results with those obtained within differ- 
ent other approaches. At moderate couplings, the 2P1 method 
gives reliable results irrespective of the temperature or the bias 
voltage both in the particle-hole symmetric and the mixed va- 
lence regimes. 

Considering the effective coupling one finds that the resum- 
mation can be thought of as using a frequency dependent four 
point vertex. One way to go beyond the approximations of 
this study is to take into account a higher loop expansion of 
the 2PI part T2{D'\ of the Hubbard-Stratonovich transformed 
theory. While this requires more numerical effort, since the 
self energies in higher loop approximation include at least one 
inner vertex, for which an integration has to be carried out, it 
is nevertheless possible with current computers. 

An expansion of the presented scheme to the non- 
equilibrium RG approach put forward in Refs. lfl6[[l7ll is un- 
der development. It was shown that in the i-channel trun- 
cation this RG scheme is equivalent to the 2P1 scheme used 
in this work. Beyond this one can take into account resum- 
mations also in the t and u channels, also called spin-singlet 
or electron-hole and electron-electron channels, respectively, 
which are expected to yield the Kondo-resonance behaviour 
for large U U^^. 
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Appendix A: Equal-time propagator 

The one-loop contribution to the self-energy gives the 
mean-field approximation ( l20l i of the dynamic equations ( fT6] |. 
Note that Dcr{t, t') has a jump at equal time arguments. Hence 
and one has the choice between D-o-{t, f + e) and D-a-{t, t - e), 
or some linear combination of them. 



Um D^t, f - e) = dAt)dl(t) = 1 - «^(f), 
lim DAt, f -H e) = -dl(t)dAt) = -n^(0 (Al) 

e->+0 

It is well known that this ambiguity in the 2 point correlation 
function corresponds to different orderings in the operator lan- 
guage. For example, choosing the second possibility in ( lAlb 



corresponds to the operator Hamiltonian 

H^Yj Eo<^4dcr + Ud\d^dldi. 



(A2) 



This can be shown by using the fact that a symmetrically or- 
dered Hamiltonian corresponds to a symmetrised path integral 
prescription. Thus, 

t U t t t t 

^/symm = 2_^Ei)i^dl.d,y + —{d^d^d ^d^ - d^^d^d'^di 

= Y^(EQ^-^Yld'T + Ud\d^dldl. (A3) 



is related to 



symm,o" 



Hm D-a^it, t- e)+ lim D_o.(f, t + e) 

e->+0 £->+{) 



F-At,t). 



(A4) 



By analyzing the equations of motion, one finds that this is 
equivalent to using the original Hamiltonian (IA2I) . and the fol- 
lowing prescription 



= -iUyi-a-6{t - t'). 



(A5) 



For higher-order contributions no such complication of order- 
ing is present. 



Appendix B: Numerical Implementation 

The naive (Euler) discretisation of the equation ( IBll l breaks 
down quickly, because the absolute value of F{t, t') is not con- 
served, even for zero right-hand side R. 



{id,-Eo)F{t,t')=R 



(Bl) 



One gets a much better discretisation if one first inserts an 
ansatz which describes the free behaviour, and then discretises 
the equation for the remaining time evolution. 

For this one introduces the new functions H{t, t') and s{t, t') 



F(t, t') = g-'^o^'-' ^H{t, t'), pit, t') = >s(t, t'). (B2) 

Without interactions, the solution is H{t, t') - const and 
s{t, t') - const. Writing the discretised equations for H{t, f), 
and similarly for s{t, t'): 

H{t + Af, t') = -ie''^°^'-''^RH{t, t')At + H(t, f), (B3) 

and multiplying this equation with e-'£o('+A/-r') ggj.^ 
equation for F{t, f): 



F{t + At, t') = -ie-'^'^'RAt + Fit, t') 



(B4) 
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This equation defines a discretisation of (IB II ) that is dif- 
ferent from the usual Euler discretisation, but has the same 
continuum limit. We see that using this algorithm, the free 
evolution is reproduced without numerical errors. Hence for a 
moderate right-hand side one can use a big timestep, and still 
conserve the unitarity of the equation to a good degree. 

Solving the 2P1 equations, we set p(f, t) = ; explicitly, as 
prescribed by the the anticommutation relations. Everything 
else is calculated using the discretisation ( IB4b of the EOM. 
On the right-hand side of the EOM, the memory integrals for 
a particular time t require field values in the past, which have 
been computed before. Therefore, our scheme allows an ex- 
plicit calculation of the next timestep from known quantities. 
We typically use Af - 0.01 - 0.02, which gives the continuum 
limit to a good accuracy. 

For zero magnetic field, when Eqi - Eqi, one can use 
a symmetric initial condition to save numerical resources. 
These symmetries apply in the second order truncation as well 
as in the resummed truncation. We discuss the resummed 
equations below: starting from an initial condition, where the 
following equations are satisfied for f = f' = 0, then they hold 
for any t, t' > 0. 



Goiit, t') = Gi2(f,' t), G„(f, f) = G22(f, t'), (B5) 
nii(f, t') = n22(f, t'), D^(t, t') = Di(f, f). 

The real constraint set by these equations is only D|(0, 0) - 
Z)|(0, 0), the rest follows automatically via the constraint 
equations dSTT i of the scalar fields and the EOM for the 
fermionic fields. 



Appendix C: Decomposition of the propagator G 

According to 0l1 l. G{t, t') is decomposed as G(f, t') - 
G(t, t') + iUcri6{t - t'), where crj is the first Pauli matrix. The 
equation for G reads 



G21 ?22 j _ ^.jy^j-j ^ ^ _ jy2 I ni2 Till 

Gil G12 / \ 1122 n2i 



(CI) 



The fermionic and bosonic self energies can be decomposed 
in the usual way into F and p parts 



F _ 1 _ 

2^0- = -^(Fa-Gfa-cT - -Pa-Gpa-a-), 



^o- - ^ 1 (Po-Gf (j-o- + Fa-Gpa-cr), 

= I77^|2_i|p^|2, 

= 2Re(0-), 



(C2) 



where, similar to the fermions, the scalar propagator is de- 
composed into 

Gijit, t') = Gf/f, f) - ^signc(f - t')GP.(t, t'). (C3) 

In this setup, the ordering problem also appears, in the equa- 
tions for the field average of the;^' fields. It can be dealt with 



similarly as in the previous subsection: symmetric Hamilto- 
nian corresponds to a symmetrised path integral. From this 
one deduces the prescription for the non-symmetrised Ander- 
son Hamiltonian ([T). Eq. ( fT6b corresponds to the ordering of 
^ when using 



(C4) 



The scalar constraint equation jClb can be decomposed using 
the identity 

(/A*B)(f,f') = r dzAP(t,z)B''{z,t') 
Jo 

dzA^{t,z)BP{z, t') 



Jo 



-signc(f 



AP{t,z)BP(zj'). (C5) 



Appendix D: Effective coupling for 5-cliannel resummation 

In this appendix we give the details of the calculation of 
the effective coupling given in Section IVIBI In the zero- 
magnetic-field case, where Eq-^ - Eqi, and the symmetries 
of (|B51 l hold, one can define: 



U 

n(f,f') = f/nii(f,f'). 



u 



(Dl) 



Hence, the constraint equations dSTT l can be written as (sup- 
pressing time arguments) 



A = /n*B-n, B = /n*A. 



(D2) 



After the transient processes are finished, in the stationary 
state, all two-point functions will only depend on the differ- 
ences of the time coordinates. In this case we can push the 
initial time to negative infinity, and we can use the following 
decomposition identities for convolution on the Schwinger- 
Keldysh contour starting at fo = 

i(X*Y)'' = X'^ * Y'" - X'" * Y"^ , 
i(X*Y)P = X'^ *YP -XP ^Y"^, 



i{X * Yf = X* * y*, 
i(X * Y)"^ = -X"^ * Y^ 



(D3) 



where the retarded and advanced two-point functions are de- 
fined as 



G^(f,f') = e(t-t')GP{t,t'), 

G^{t,t') = 9(t' -t)GP{t,t'). 



(D4) 



Using this in Eq. ( ID2l i. after eliminating B and Fourier 
transformation with respect to f - f', we obtain 

A^ = n^n^A'" - n^n'^A'^ -Yf + n^'n^'A'*, 

AP = n^n^A'^ - ifliPA^ -UP + TiPlfA^ (D5) 
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where we have omitted arguments (oj). With the help of ( ID3I ). where 
Eq. ( ID2I ) can be rewritten as 

(i + n«A^-n^A^)(n^n^-i) = n«n^-i. (D6) , , i-n^n^ i - P^l- 

^ ' ^e.t<'-) = ((n«)2-i)((m)2-i)^(n«)2-ip- ^^^^ 

We can write 

A'' - -n^'f/eff. a'' - -YfUf^s (D7) which is the expression given in Eq. dSTl ). 
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